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We report a measurement of the process 77* — > tt° with a 759 fb _1 data sample recorded with 
the Belle detector at the KEKB asymmetric-energy e + e~ collider. The pion transition form factor, 
F(Q 2 ), is measured for the kinematical region 4 GeV 2 ~ Q 2 ~ 40 GeV 2 , where -Q 2 is the invariant 
mass squared of a virtual photon. The measured values of Q 2 \F(Q 2 )\ agree well with the previous 
measurements below Q 2 ~ 9 GeV 2 but do not exhibit the rapid growth in the higher Q 2 region 
seen in another recent measurement, which exceeds the asymptotic QCD expectation by as much 
as 50%. 
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I. INTRODUCTION 



Recently, the BaBar collaboration reported a mea- 
surement of the photon-pseudoscalar meson (P) 
transition form factors (TFF) from the process 
77* — > P in the reaction e + e~ — ¥ (e)eP, where 
P = 7T° [T] (Fig. [TJ, 7/ or 7/ [2], and (e) and e de- 
note an untagged and tagged electron/positron, re- 
spectively. The momentum transfer Q 2 (where —Q 2 
is the squared invariant mass of a space-like photon) 
explored by BaBar extends up to ~ 40 GeV 2 com- 
pared to the previous measurements by CELLO for 
Q 2 < 2.5 GeV 2 [3] and by CLEO for Q 2 < 9 GeV 2 g]. 
While the Q 2 dependence of TFFs measured by 



BaBar for the 77 and 77' agrees fairly well with theoret- 
ical expectations, that for the ir° shows rapid growth 
for Q 2 > 10 GeV 2 . 



Lepage and Brodsky pointed out the importance 
of TFF measurements more than 30 years ago and 
proposed tests of perturbative QCD (pQCD) in var- 
ious exclusive processes [5]. We focus our discussion 
on the TFF for the pion; that for the 77 or 77' is quite 
similar. To leading order in the QCD coupling con- 
stant, a s (Q 2 ), the TFF can be written in the factor- 
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FIG. 1: A Feynman diagram for 77* — > qq — > ir in 
e + e~ collisions 

ized form [6] 

F(Q 2 ) = f dxT H (x,Q 2 ,^,a s ( f i)) 

Jo 

X</) n (x,fi), (1) 

where f n ~ 0.131 GeV is the pion decay con- 
stant, x is the fraction of momentum carried by 
a quark (q) or antiquark (q) in the parent pion, 
Th(x, Q 2 1 A*) a s(/-i)) is the hard-scattering amplitude 
for 77* — > qq (Fig.[lJ, [i is the renormalization scale, 
and <j)„(x,fj,) is the leading-twist pion distribution 
amplitude (DA) at x and fi. Note that the pion 
DA is a universal function and, once determined, de- 
scribes leading terms in all hard exclusive processes 
involving a pion. The most prominent examples are 
the pion electromagnetic form factor and the weak 
semileptonic decay rate, B ± — > ifii^vt [BJ. 

At sufficiently high Q 2 , Eq. can be written 
as [B] 

F(Q 2 ) = ^f Q dx^- + O (1/Q 4 ) . (2) 

Because the asymptotic pion DA is <jf^ y = 6x(l — 
x), it follows that Q 2 F(Q 2 ) at Q 2 — >• 00 is exactly 
normalized as [5] 

Q 2 F(Q 2 ) = V2U - 0.185 GeV. (3) 

The measurements reported by CELLO [3] and 
CLEO [4 support this prediction; the measured val- 
ues seem to approach this asymptotic value from 
below. The recent measurement by BaBar [T] has 
changed this situation drastically. The measured 
Q 2 \F(Q 2 )\ values seem to increase with Q 2 , exceed- 
ing the QCD asymptotic value by as much as 50%; 
this has created excitement and renewed interest in 
this quantity. 



Many theoretical constructs have been made to 
explain this Q 2 dependence. Early attempts pointed 
out that the BaBar result requires a nearly flat pion 
DA (e.g., jTJ [5]), which implies that the pion be- 
haves as a point-like particle. The conclusions of the 
theoretical analyses can be divided into two groups: 
those that can explain the BaBar result (e.g., [6j|9]) 
and those that cannot (e.g., |10l \TT\ ). 

This theoretical controversy seems to have been 
settled by a recent analysis of Brodsky, Cao and de 
Teramond [12 , who investigated TFFs of the it , 
T) and rj using four typical models of the meson 
DA. They concluded that the rapid growth of the 
7T° TFF with Q 2 at BaBar, together with the gentle 
Q 2 -dependence of the 77 and f/ TFFs also measured 
at BaBar are difficult to explain within the present 
framework of QCD. They also pointed out drawbacks 
of earlier theoretical attempts that explained the Q 2 
dependence of the ir° TFF at BaBar and stated that 
the BaBar tt° result, if confirmed, would imply new 
physics beyond standard QCD. 

In this Article, we report an independent measure- 
ment of the tt° TFF in the Q 2 region of interest. 



II. EXPERIMENTAL APPARATUS AND 
EVENT SELECTION 

We briefly describe the Belle detector and then 
give a description of the event selection. Finally, 
we present some event distributions to illustrate the 
selection criteria. 



A. The Belle detector and data sample 

We use a 759 fb _1 data sample recorded with 
the Belle detector at the KEKB asymmetric-energy 
e + e~ collider [L3J . We combine data samples col- 
lected at several beam energies: at the T(45) res- 
onance (y/s = 10.58 GeV) and 60 MeV below it 
(637 fb _1 in total), at the T(3S*) resonance (y/s — 
10.36 GeV, 3.2 fb _1 ) and near the T(55) resonance 
(Vs = 10.88 GeV, 119 ftr 1 ). When combining the 
data, the slight dependence of the two-photon cross 
sections on beam energy is properly taken into ac- 
count. 

This analysis is performed in the "single-tag" 
mode, where either the recoil electron or positron 
(hereafter referred to as electron) alone is detected. 
We restrict the virtuality (Q 2 ) of the untagged-side 
photon to be small by imposing a strict transverse- 
momentum balance with respect to the beam axis 
for the tagged electron and the final-state neutral 
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pion. In this Article, we refer to events tagged by 
an e + or an e~ as "p-tag" (positron-tag) or "e-tag" 
(electron-tag) , respectively. 



1. Belle detector 

A comprehensive description of the Belle detector 
is given elsewhere |14| . We mention here only the 
detector components essential for this measurement. 
Charged tracks are reconstructed from the drift-time 
information in a central drift chamber (CDC) lo- 
cated in a uniform 1.5 T solenoidal magnetic field. 
The z axis of the detector and the solenoid are along 
the positron beam direction, with the positron beam 
pointing in the — z direction. The CDC measures 
the longitudinal and transverse momentum compo- 
nents, i.e., along the z axis and in the rip plane per- 
pendicular to the beam, respectively. Track trajec- 
tory coordinates near the collision point are provided 
by a silicon vertex detector (SVD). Photon detec- 
tion and energy measurements are performed with a 
CsI(Tl) electromagnetic calorimeter (ECL). Electron 
identification (ID) is based on E/p, the ratio of the 
calorimeter energy and the track momentum. We 
employ this rather simple variable because the back- 
ground from hadrons and muons is small compared 
to the dominant radiative Bhabha and two-photon 
background contributions. 



2. Data sample 

To be recorded, events must satisfy one of the 
two mutually exclusive triggers based on ECL en- 
ergy: the HiE (High-energy threshold) trigger and 
the CsiBB (ECL-Bhabha) trigger [15]. The HiE trig- 
ger requires that the sum of the energies measured 
by the ECL in an event exceed 1.15 GeV but that the 
event not be Bhabha- like; the HiE trigger is vetoed 
by the CsiBB trigger. The CsiBB trigger is designed 
to identify back-to-back Bhabha events |TS] and is 
prescaled to reduce their high rate. For the purpose 
of monitoring trigger performance, the events trig- 
gered by the CsiBB are recorded once per 50 CsiBB 
triggers (i.e., are prescaled by a factor of 50). This 
prescaled event sample has been extensively used in 
our analysis to select a subsample of signal events, 
and to calibrate and test the trigger simulation for 
the CsiBB veto implemented in the HiE trigger. 

We do not use information on track triggers be- 
cause they require two or more charged tracks; our 
signal has only one charged track. 



B. Single-track skim 

To simplify the logistics of the data analysis, we 
use a "single-track" skim whose requirements are the 
following. 

(51) The event is triggered by the HiE or the CsiBB 
trigger. 

(52) There is only one track that satisfies p t > 
0.5 GeV/c, -0.8660 < cos (9 < 0.9563, dr < 
1 cm and \dz\ < 5 cm. There are no other 
tracks that satisfy p t > 0.1 GeV/c and dr < 
5 cm, \dz\ < 5 cm in the entire CDC volume. 
Here, p t is the transverse momentum in the 
laboratory frame with respect to the positron 
beam axis, 9 is the polar angle of the momen- 
tum direction with respect to the z axis, and 
(dr, dz) are the cylindrical coordinates of the 
closest approach of the track to the beam axis. 

(53) There are one or more neutral clusters in 
the ECL, whose energy sum is greater than 
0.5 GeV. 

These conditions are efficient in selecting a signal 
process within the kinematical regions of e + e _ — > 
(e)e7r° in which one electron escapes detection at 
small forward angles. They also select radiative 
Bhabha (e + e~ — > (e)ej) events with a Virtual 
Compton (VC) process configuration, in which one 
of the electrons escapes detection in the forward di- 
rection; such events are used for trigger calibration. 

C. Selection criteria for signal candidate events 

A signal event consists of an energetic electron and 
two photons as described below. Candidate events 
are selected from the single-track skim with the fol- 
lowing selection criteria. The kinematical variables 
are calculated in the laboratory system unless other- 
wise noted; those in the e + e~ center-of-mass (cm.) 
frame are identified with an asterisk. 

(1) Events are triggered by the HiE trigger. We 
also use events triggered by the CsiBB trigger, 
separately, for the Q 2 = 4 — 6 GeV 2 region of 
the e-tag. 

(2) There is exactly one track as required in the 
skim condition (S2). 

(3) For electron ID, we require E/p > 0.8 for the 
candidate electron track. 

(4) The absolute value of the momentum of the 
electron is greater than 1.0 GeV/c, where the 
electron energy is corrected for photon radia- 
tion or bremsstrahlung in the following way. In 
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a 3° cone around the track, we collect all pho- 
tons in the range 0.1 GeV< E 1 < p e /3, where 
p e is the measured absolute momentum of the 
electron track. The absolute momentum of the 
electron is replaced by p e + HE 1 . 

(5) We require exactly two photons above 0.2 GeV, 
excluding those used for bremsstrahlung and 
final state radiation recovery in (4). We des- 
ignate the energies of the two photons by 

for the higher and E~ t 2 for the lower. 

(6) When there are one or more extra photons with 
energy greater than 0.05 GeV, we require that 
there be no ir° candidate reconstructed with a 
two-photon combination different from the pair 
in (5) in the following way. We search for a 7r° 
candidate reconstructed from two photons that 
satisfies a goodness-of-fit criterion, \ 2 < 9, in 
a 7r°-mass-constrained fit. If we find such a ir° 
(including the case when only one photon in 
(5) is shared in the 7T°), the event is rejected. 

(7) The two-photon energy sum i? 77 defined as 

-E 77 = Ery\ + Ej2 (4) 

satisfies E 11 > 1.0 GeV. 

(8) The cosines of the polar angles for the electron 
(6 e ) and the two photons ((9 7 i and 6^2) must 
be within the range -0.6235 < cos 6^ < 0.9481 
(i = e, 7I or 72). This is the sensitive region 
for the HiE trigger. 

(9) We reject events with a polar-angle combina- 
tion (cos6* e , cos# 77 ) that falls in the predeter- 
mined two-dimensional angular area shown in 
Fig. [2] where 6* 77 is the polar angle of the mo- 
mentum summed for the two photons. This 
rejection is referred to as the "Bhabha-mask" 
rejection. The Bhabha mask is introduced to 
limit the fiducial region to a range where the ef- 
ficiency of the HiE trigger is well determined. 
The details of the Bhabha mask are given in 
Appendix [X] 

(10) The energy asymmetry for the two photons, 
defined as 

_ E yl - E, (2 

-^asym — ^ > \<J ) 

satisfies -E a sym < 0.8 to eliminate low-energy 
background photons. 

(11) The polar angle difference, AO — \0~i — 9^1 for 
the two photons multiplied by i? 77 satisfies 

E 17 A9 > 0.18 radian • GeV. (6) 

This rejects the large background from ra- 
diative Bhabha events with a conversion, i.e. 



e + e~ — > (e)e7, 7 — > e + e~ , where the electron- 
positron pair from a photon conversion tends 
to separate in a common plane perpendicu- 
lar to the magnetic field (which is along the 
z axis), and the e + e~ pair is misreconstructed 
as a 77 candidate. 

(12) We exploit the charge-direction correlation 
("right-sign") for the single-tag two-photon pro- 
cess between the electric charge of the tagged 
electron (<?tag) and the z component of the mo- 
mentum of the e77 system in the e + e~ cm. 
frame by requiring 

<?ta S x (p;, e +p; iry ) <o. (7) 

(13) We require 0.85 < £ ra tio < 1.1, where 

^*measurcd 
^ratio = ^expected ' ( 8 ) 

This requirement is motivated by a three-body 
kinematical calculation for e + e~ — > (e)e7r°, im- 
posing four-momentum conservation condition 
Pinitiai(e + e _ ) = Pfinai((e)e77) wherein the di- 
rection of the 7T° momentum is taken to be par- 
allel to that of the observed 77 system in the 
e+e~ cm. frame. The expected E* value is 
obtained by assigning the nominal ir° mass to 
the 77 system. This calculation is not sensitive 
to the input 77 invariant mass. 

(14) We reject events with a back-to-back configura- 
tion in the e + e~ cm. frame, to eliminate back- 
grounds from Bhabha events in which a track is 
not reconstructed; we require C*( e :77) < 177°, 
where C*( e >77) i s th e opening angle between 
the electron and the 77 system. 

(15) We require the electron and the 77 system 
to be nearly back-to-back when their momen- 
tum directions are projected onto the r<p plane: 
a(e,77) < 0.1 radians, where the acoplanarity 
angle a is defined as the difference between the 
opening angle of the two directions on the pro- 
jected plane and ir. 

(16) We require transverse momentum balance in 
the e+e~ cm. frame, |Sp*| < 0.2 GeV/c, 
where 

|5X| = |p* e +p* 77 |. (9) 

With the above selection criteria, we collect 
e + e~ —¥ (e)e77 candidates in the broad 77 invariant- 
mass range of 0.07 GeV/c 2 ~ M 77 ~ 0.50 GeV/c 2 . 
The minimum 77 invariant mass is determined 
mainly by selection criteria (7) and (11) on £? 77 and 
E-yjA0. We do not apply any explicit selection on 
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FIG. 2: Two-dimensional angular regions in (cos^e, cos# 77 ) defined as the Bhabha mask (see Appendix |A| for 
the details). Events in the yellow (white) regions are selected (rejected). The gray region is outside of the angular 
acceptance. The mask is symmetric for e and 77. The scattered dots in the figure show signal Monte Carlo simulation 
events without the Bhabha veto applied. 



M 77 for the determination of the ir° signal yield. In- 
stead, we later obtain the 7r° yield by fitting the M 77 
distribution in Q 2 bins separately for the p-tag and 
e-tag samples: the selection requirements for the two 
tags are very different because the beam energies and 
the detector are asymmetric. 



D. Event distributions 

In this subsection, we illustrate how the signal 
candidates are selected. Figure [3] shows the energy 
asymmetry (i^asym; defined in criterion (10)) of the 
two photons for events in which the electron and 
photons satisfy the selection criteria for momenta, 
energies and angles. Events with a large asymme- 
try are suppressed by the minimum photon energy 
requirement (0.2 GeV), yet a peak at E asym > 0.8 
is still visible in the data from low-energy photon 
background (Fig. 3^a)). 

Figures [4] and [5 show -E ra tio versus M 77 and |Ept|, 
respectively, for the samples after applying selection 
criteria (11) and (15) for E^AO and the acopla- 
narity angle, respectively. The latter is plotted 
only for events in the vicinity of the pion mass, 



0.120 GeV/c 2 < M 77 < 0.145 GeV/c 2 , where the 
-Eratio requirement is not applied. We find a sig- 
nificant enhancement in the signal region of the 
experimental data near -E ra tio=l, M T7 = m, and 
|Spjj| = 0, where is the nominal tt° mass. A nar- 
row enhancement in the experimental data around 
-EVatio = 1-15 (Fig. |4ja)) is due to background from 
radiative Bhabha events and is reproduced by the 
background Monte Carlo (MC). The signal in the 
experimental data is well separated from the low- 
energy 7T° background in the region near | Ep^ | = 
(Fig. [5}. 



In Fig. [6] we show the scatter plot for E 17 A9 vs. 
M~ n for the final samples in the range 15 GeV 2 < 
Q 2 < 40 GeV 2 , where only the corresponding selec- 
tion criterion is not applied. The large background 
from radiative Bhabha events is present primarily in 

the range 0.05 GeV/c 2 ~ M 77 ~ 0.20 GeV/c 2 below 
the selected boundary in E^AO, with no contami- 
nation seen above the boundary. This background is 
discussed in more detail in Sec. IIVBI 
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FIG. 3: Energy asymmetry (-E asym , denned in selection criterion (10)) for the two photon candidates for events in 
which p e , _E 77 and # 77 satisfy the criteria for (a) experimental data and (b) signal Monte Carlo samples. The dashed 
lines show the selection boundaries. Arrows indicate the accepted regions. 
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FIG. 4: Scatter plots for -Eratio vs. M 77 for (a) experimental data and (b) signal MC samples. See the text for the 
selection conditions. The horizontal band between the dashed lines in (b) is the accepted region. 



III. SIGNAL MONTE CARLO AND 
ESTIMATION OF THE DETECTION 
EFFICIENCIES 



In this analysis, MC programs play an essential 
role in defining the selection criteria and determining 
the detection efficiencies. In this section, a specially 
developed signal MC program is described. The ef- 
ficiency determination is then described. A detailed 
discussion of the trigger simulator tuning is given 
in Appendix [B] where radiative Bhabha events are 
compared with MC events generated by Rabhat [TB] . 
This MC is a specially developed program to deal 
with the VC process in which an electron is scat- 
tered forward and remains undetected. 



A. Signal Monte Carlo, TREPSBST 

To calculate the conversion factor between the 
transition form factor and the differential cross sec- 
tion, to estimate the overall efficiency using signal 
MC events and to perform other studies, an MC 
program, TREPSBST, was developed for the signal 
process e + e~ — > (e)eir°. 

TREPSBST implements the formulae of Eqs. (2.1) 
and (4.5) in Ref. [T7] and is based on MC pro- 
gram TREPS [IB], which was modified to match the 
condition that the equivalent photon approximation 
(EPA) is not applied. In TREPSBST, we always 
retain the condition Q\ > Q\ for virtuality of the 
two colliding photons by requiring Q\ > 3.0 GeV 2 
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FIG. 5: Scatter plots for -B ra tio vs. |Ep*| for the (a) experimental data and (b) signal MC samples. See the text for 
the selection requirements. The solid box in (b) shows the selection region. 
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FIG. 6: Scatter plots for E-,-,A8 vs. M 77 for (a) the experimental data and (b) signal MC samples. This distribution 
is for the final samples in the range 15 GeV 2 < Q 2 < 40 GeV 2 where only the _E 77 A# criterion is not applied. The 
dashed lines show the selection boundaries. Arrows show the accepted regions. 



and Q 2 < 1.0 GeV 2 for the ranges of integration and 
event generation. 

This Q\ range is sufficient to generate signal events 
over the kinematical region for the selection 
criterion, and is also used to define the differential 
cross section da/dQ 2 for the signal process in this 
analysis (see Sec. VC). 



The program calculates the differential cross sec- 
tion da jdQ\ for the given Q\ points and, separately, 
generates MC events with a specified Q\ distribu- 
tion. The shape of the distribution, in principle, does 
not affect the efficiency determination, which is eval- 
uated for each Q\ bin with a bin width of 1 GeV 2 . 
However, we assign a systematic error from the un- 
certainty in the Q 2 -dependence effect for some bins 



in the lowest Q 2 region (see Sec. VEl. 

The differential cross section for the signal process 
is calculated at tree level assuming a g 2 coupling pa- 
rameter derived from T 77 = 8 eV (defined for real 
photons as in Eq. ( |T2| ) below). Both the transition 
form factor squared — fi 7, (Q 2 )| 2 and differential cross 
section da/dQ 2 are proportional to the parameter 
9 2 



2 and thus their ratio 



2A(Q 2 ) 



da/dQ 2 

WWW 



(10) 



which has to be used to convert the measured differ- 
ential cross section to the form factor, is not affected 
by the assumed T 77 value [T7] . In Eq. ( 10 1 , we explic- 



itly include a factor of two to show the contributions 
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from both p- and e-tags to the coefficient. The tran- 
sition form factor F(Q 2 ) is defined as 



\F{Q 



2\|2 



lim |F(Q 2 ,Q 2 



2\|2 



(11) 



The normalization of the form factor F(Q 2 : Q%) is 
defined based on the relation in Ref. H], where it is 
extrapolated to the real two-photon coupling (77 — » 



|F(0,0)| 2 = 



1 647Tl\ 



(Ana) 2 



(12) 



here a is the fine structure constant and T 77 is the 
two-photon decay width for the neutral pion. 

A calculation of the conversion factor 2A(Q 2 ) re- 
quires a 3-fold numerical integration, which is per- 
formed over the following variables: Q 2 , the photon 
energy of the small-Q 2 side (u>2), an d the acopla- 
narity for the two colliding virtual photons (A(p). 
The choice of kinematical variables is discussed in 
Ref. PI]. 

We use the following assumption to extend the 
form factor to non-zero Q 2 : 



\F(QlQl 



TO?, Q)l 
i + S 



(13) 



where m p = 0.77 GeV/c 2 . The ratio 2A(Q 2 ) is cal- 
culated for the Q 2 range of the present measurement 
with a step size of 1 GeV 2 . Convergence and stabil- 
ity of the numerical integration have been carefully 
verified to an accuracy of ~ 10~ 3 . We have applied 
spline interpolation for 2A(Q 2 ) when necessary. The 
results obtained for 2A(Q 2 ) are plotted in Fig. [7] for 
s/s = 10.58 GeV, Q^ max = 1.0 GeV 2 and to leading 
order in QED. 

The event generation is performed using the same 
integrand but including initial-state radiation (ISR) 
effects from the tag-side electron. Inclusion of ISR 
changes the kinematics and Q 2 significantly. Mean- 
while, ISR from the untagged side has little ef- 
fect because ISR is nearly parallel not only to the 
initial-state electron but also to the final-state un- 
tagged electron. We use an exponentiation tech- 
nique (20) for the photon emission based on the pa- 
rameter i] = (2a/7r)(log(Q 2 /TO 2 ) — 1) and the prob- 
ability density for the photon energy distribution, 
dP(rk) oc r^T 1 drk, where 



771* 

= ^ISR 

" E* 



(14) 



As an approximation, the photon is always emitted 
along the electron direction. We set the maximum 




30 

q2 (QeV2) 

FIG. 7: Calculated 2A(Q 2 ) function for ^/s = 10.58 GeV. 



fractional energy of radiation r™ ax = 0.25. Because 
events with a large rj,, typically above 0.1, are re- 
jected by selection criterion (13) using -E^tio (see 
Sec. ~ 



VE2 



and Fig. 23 1 



MC 

generation is large enough to generate MC events 
over the selected range. In this configuration, the 
correction factor 1 + 5 to the tree-level cross section 
from the untagged side is close to unity [3T]. Radia- 
tion of low-Q 2 virtual photons is strongly suppressed 
by the emission of ISR in the untagged side. 

We generate events with a virtuality of the tagged- 
side photon Q\ distributed with a constant form fac- 
tor. The Q 2 of each event is modified by ISR. We 
use the momentum of the ISR photon to estimate 
the true Q 2 value in signal MC events and to study 
the Q 2 dependence of the detection efficiency. 



B. Efficiency determination 

The selection efficiency is calculated using the sig- 
nal MC, the detector simulation (GSIM) and the 
trigger simulation (TSIM). The efficiency is esti- 
mated for each Q bin for p- and e-tags separately. 
We show the Q 2 dependence of the efficiency for a 
typical run at the T(4S") energy and another at the 
T(5S*) energy in Figs. [8] and [91 The overall detec- 
tion efficiency shown in Fig. [8] includes the trigger 
efficiency and the Bhabha-mask selection. The trig- 
ger efficiency (Fig. [9]) is shown within the acceptance 
region after the Bhabha-mask selection is applied. 

A measurement corresponding to these efficien- 
cies is reported only in the region 6 GeV 2 < Q 2 < 
40 GeV 2 , where the HiE trigger with the Bhabha 
mask is used. The complicated shapes are due to 
the structure of the Bhabha mask. The efficiencies 
are obtained for each Q 2 bin, where the value of Q 2 is 
evaluated using generated four-momentum vectors of 
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signal MC events with ISR taken into account. We 
generated MC samples of 8 x 10 6 events at T(45) 
and 2 x 10 6 events at T(5S); these samples corre- 
spond, approximately, to an integrated luminosity a 
few dozen times larger than our data. The efficiency 
obtained in bins of Q 2 is summarized in Table |lj 

To validate the efficiency calculation, we compare 
the distributions for the 7r° laboratory angle in sig- 
nal candidate events with those from the sign al MC 



Figure 11 shows scatter plots of Q 2 vs. M 77 for 



samples in different Q regions in Fig. 10 Data 
events are shown in the range of two-photon invari- 
ant mass 0.120 GeV/c 2 < A/ 77 < 0.145 GeV/c 2 af- 
ter subtracting the background that was estimated 
using the 7r°-mass sideband region 0.170 GeV/c 2 < 
M 77 < 0.195 GeV/c 2 using the fitting procedure de- 
scribed in Sec. |IV A} Qualitative agreement between 
the data and MC in Fig. [10] confirms the validity 
of the efficiency calculation based on the signal MC. 
The validity of this calculation has also been checked 
in data using a sample of radiative-Bhabha events, 
which have no trigger bias. The details are described 
in Appendix IB] 



IV. EXTRACTION OF THE tt° YIELD AND 
BACKGROUND ESTIMATION 

In deriving the differential cross section, da/dQ 2 , 
we first extract the 7r° yield in each Q 2 bin by fitting 
the M 77 distribution. Possible backgrounds are then 
identified, evaluated and subtracted. These analysis 
steps are described below. 



A. Extraction of the 7r° yield 

The 7T° yield is extracted in each Q 2 bin for the p- 
tag and e-tag samples. For each signal candidate 
event, Q 2 is calculated using the measured four- 
momentum of the detected electron (p e ) by 

Qrcc = -Ob earn Pe) 

= 2E^ mn E* e (l + q tag cose* e ), 



(15) 

where Pbcam is the nominal four-momentum of the 
beam particle with the same charge as the detected 
electron; the right-hand side is given by the beam 
energy E! eam and the observables of the tagged elec- 
tron in the e + e~ cm. frame. We do not apply a cor- 
rection for ISR on an event-by-event basis; instead, 
this effect is taken into account in the Q 2 unfolding 
described in Sec. IV Al 



the final experimental samples. After dividing them 
into Q 2 bins, we fit the M 77 distribution in the region 
0.08 GeV/c 2 < Af 77 < 0.30 GeV/c 2 to determine the 
7T° yield. 

To extract this yield, we fit the two-photon invari- 
ant mass distribution with a double Gaussian de- 
scribing the signal and a quadratic polynomial ap- 
proximating the background, namely 



\-bx + 
1-r 



A 



2-ku 



2(kcr) 2 



{re 

2 

1, 



(16) 



where x = Af 77 , a, b and c are the background pa- 
rameters, A is signal yield, and m and er are the mass 
and width (i.e., mass resolution), respectively, for 
the signal. For Q 2 bins with poor statistics, the last 
two fit parameters are fixed to values estimated from 
signal MC, as described below. The parameters r, k 
and Am approximate the relative size, width and po- 
sition of the two Gaussian components, respectively, 
and are fixed to values from MC simulation in a Q 2 - 
dependent manner. 

In Q 2 bins with low signal statistics, the determi- 
nation of the width is poor and is significantly corre- 
lated with the yield. For such bins, we fix the mass 
and width parameters to values that are expected 
from the signal MC, taking into account a small de- 
viation measured using the data: the difference be- 
tween data and MC is 6m ~ — 1 MeV/c 2 , whereas 
the mass resolutions show good agreement. These 
parameters are fixed in the fit for Q 2 > 20 GeV 2 for 
the p-tag, and Q 2 = 7 - 8 GeV 2 or Q 2 > 12 GeV 2 
for the HiE and Q 2 < 6 GeV 2 for the CsiBB for 
the e-tag. The systematic uncertainty induced by 
fixing these parameters is dominated by the large 
statistical error on the yield difference between the 
constrained and non-constrained fits. We estimate 
this uncertainty by comparing the results obtained 
using two different sets of fit functions, as described 
in Sec. IVEl 

The following values for the fixed parameters were 
obtained from the MC: m = 134 - 137 MeV/c 2 , 
a = 6 - 9 MeV/c 2 , Am = -5.5 - 0.0 MeV/c 2 , 
1 — r = 0.11 — 0.42 depending on Q 2 and the p- 
or e-tag; k is fixed at 2.4. The absolute values for 
these parameters tend to be larger when Q 2 is large. 
There are no significant anomalous tails in the sig- 
nal M 77 distribution down to the 0(1O -2 ) level of 
the peak, according to the signal MC, as shown in 
Fig. [T2} Figure [13] shows the results of the fit to data 
in five representative Q 2 bins. The yields of the fits 
are summarized in Table |T] and are labeled N rec . 
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FIG. 8: The total selection efficiencies as a function of Q 2 determined by MC for runs at the T(4S) (solid histogram, 
where the error bars show statistical errors) and T(5S) energies (dashed histogram). 
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FIG. 9: The trigger efficiencies determined from MC for a typical run at the T(4S) energy. The error bars show 
statistical errors. 



1. Backgrounds from e + e annihilation processes 



B. Estimation of background contamination 

The measured 7r° yield could contain some con- 
tribution from background processes in which real 
7r°'s are present in the final state. In this subsection, 
we explain how these resonant backgrounds are esti- 
mated and subtracted. First, possible backgrounds 
from e + e~ annihilation processes are examined. A 
potentially large background from radiative Bhabha 
events is also checked. Finally, we study the back- 
ground contribution from the processes 77* — > 7r°7r° 
and e + e~ — > (e)ep /oj, p° /u — > ir°j. 



Backgrounds from e + e~ annihilation are esti- 
mated by examining the yield of wrong-sign events, 
in which the tagged electron has the opposite sign to 
that expected in selection criterion (12). The elec- 
tron charge in the final state of the annihilation pro- 
cesses should have no significant forward-backward 
asymmetry. The contamination of this kind of back- 
ground in the signal sample is expected to be the 
same as the number of events found in the wrong- 
sign samples. 

We have studied the wrong-sign events by examin- 
ing the parameter £ r atio> whose distribution is shown 
in Fig. 14 'a). We find that a small number of events 
peak near -E ra ti = 1 for the wrong-sign sample. 
However, there is no visible ir° peak in the M 77 dis- 
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TAB LE I: Numbers of reconstructed events (JV r0 c) obtained from the fit to the experimental M 11 distribution 



(Sec. IV A I, the number of events after the unfolding (N cor ), background fractions and efficienc ies for each Q bin, 
listed separately for p- and e-tags. The first column is the representative Q 2 for that bin (see Eq. ( 18 1). N Tec and N COI 
include both signal and background yields. Note that the "HiE + 50*CsiBB" samples are used for the Q 2 = 4-6 GeV 2 
bins in the e-tag case. 
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tribution in any Q 2 regions as shown in Fig. 14 d). 
The yield in the tt° peak is 1.2±0.9 events in the sum 
of p- and e-tags with Q 2 < 40 GeV 2 for the total sam- 
ple. We conclude that the neutral pions from e + e~ 
annihilation processes do not contaminate the sig- 
nal region -E ra ti — 1, and the main background near 
£-ratio = 1 is from the QED process e + e~ — > e + e~77, 
which is of order a 4 . Although this process has a 
charge asymmetry, this is not a problem for back- 
ground estimation since it does not peak at the pion 
mass. 

We also conclude that any hadron or muon 
misidentified as a tagged electron is not accompa- 
nied by neutral pions and does not contribute to 
the peaking backgrounds. Hadrons or muons would 
appear with the same probability in the wrong-sign 
sample if they were significant in the right-sign sig- 



nal sample. As seen in Fig. 14 b,c), the tail of the 
-E-ratio distribution near the lower selection boundary 
(0.85) is dominated by the ISR tail of the signal pro- 
cess, which is accounted for in the signal MC down 
to rfc = 0.25. Some contributions from the pe aking 



backgrounds ((e)e7r°X) discussed in Sec. 
4 will also be present. 



IV B 



3 and 



2. Radiative Bhabha backgrounds 

Estimation of the radiative Bhabha backgrounds 
where the photon is misreconstructed neutral 
pion is challenging. While the probability of such 
misreconstruction is small, the cross section of this 
background process is much greater than that of 
the signal process. A misidentification rate as small 
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OttO (deg) 

FIG. 10: Comparison between the n° angular distri- 
butions in signal candidate events in data (circles) and 
signal MC (histograms) for different Q 2 regions. Closed 
(open) circles and solid (dashed) histograms are from 
events of p-tag (e-tag). In the data, the background con- 
tribution is subtracted. The MC events are normalized to 
the experimental yield after the background subtraction 
in each Q 2 region. 



lution satisfying the required three-body signal kine- 
matics even if the photon is emitted at a forward 
angle. However, in this case, the background would 
have a broad -E-Vatio distribution concentrated at zero 
because the secondary ir° 's are dominantly of low en- 
ergy. There is a depletion of events around 0.8 in the 
-Bratio distribution as shown, e.g., in Fig. |14[a) , even 
for the high Q 2 regions; we conclude that the contri- 
bution from this type of background is negligible. 

There is a small probability that a photon from a 
radiative Bhabha converts to a neutral pion at the 
beam pipe or the central part of the detector, where 
the material thickness is kept to a minimum, via a 
soft nuclear interaction -fA —> ir°A'. The pion is 
produced at forward angles with a small momentum 
transfer; the nucleus A can break up with an en- 
ergy slightly larger than the binding energy. This 
cross section is poorly known, but we estimate that 
it is smaller than that for the elementary process 
7P — > 7r p per nucleon. According to the previous 
measurements [22], the cross section for scattering 
angles below 10° is about 0.04 for incident pho- 
ton energies between 2 and 9 GeV. This cross sec- 
tion corresponds to a probability of 1.0 x 10~ 7 for 
a photo-nuclear conversion to a neutral pion with a 
similar energy. From the observed number of radia- 
tive Bhabha events, we estimate that the contam- 
ination is at most 0.5 event in the highest Q 2 bin 
and negligible for all other bins. We thus conclude 
that the background from this source in all Q bins 
is small, and assign 3% and 1% uncertainties only to 
the two highest Q 2 bins. 



as C(10 -6 ) can lead to background contamination. 
However, as described below, even in a conservative 
estimate, the background is negligibly small, safely 
less than one event in each high Q 2 bin. 

The photon conversion process, 7 — > e + e~ , is the 
main source of radiative Bhabha background, but 
almost all such events are rejected by selection cri- 
terion (11). Moreover, conversions can only give a 
broad distribution in M 77 . We estimate the back- 
ground level to be negligibly small for Q 2 < 20 GeV 2 



by extrapolating the E 11 /S.6 distribution. Figure 15 



shows this distribution for the experimental data 
and for signal and background MC, where the back- 
ground contribution falls exponentially while the sig- 
nal has a nearly flat distribution. 

To check for additional background arising from 
radiative Bhabha events, we study the possibility 
that some pions from a secondary interaction or 
beam background overlap with a radiative Bhabha 
event where a photon escapes detection (e + e~ — ¥ 
(e)e(7) + 7T° (overlap)). Note that there can be a so- 



3. Background from 77* — > ir tt 

Processes with one extra particle, i.e., e + e~ — > 
(e)eir°X, where X is another ir° or photon, could 
also be a source of resonant background that would 
have to be subtracted from the measured 7r° yield in 
each bin of Q 2 . 

This type of background may have both a charge 
asymmetry that favors the right-sign and a peak at 
the 7T° mass. For large X invariant masses, this back- 
ground is suppressed by the -E ra tio selection criterion. 
However, if the invariant mass and the momentum of 
X are small, -E ra tio approaches unity, and we cannot 
use it to distinguish between signal and background 
processes. 

To study such background, we use a rather loose 
set of selection criteria to positively identify e + e~ — > 
(e)e7r°7r° candidates that contribute in a wide kine- 
matical range. There should be four or more pho- 
tons in the final state; we accept the extra pion with 
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FIG. 11: Scatter plots of Q 2 vs. M 77 for the final experimental samples, shown separately for the p-tag (left) and 
e-tag (right) events. 
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FIG. 12: M 77 distributions for the signal MC events for the selected Q 2 ranges for the p-tag (left, 7 GeV 2 < Q 2 < 
8 GeV 2 ) and e-tag (right, 20 GeV 2 < Q 2 < 25 GeV 2 ) samples. The curves are fits with a double-Gaussian function 
without background terms. 



a loose energy constraint. We apply a selection re- 
quirement on the energy of the 7r°7r° system using 
an energy ratio similar to -E ra tio to enhance three- 
body kinematical configurations. We take the ratio 
of the observed energy to the expected energy calcu- 
lated using the observed invariant mass for the 7r°7r° 
system. 

An MC generator was developed for this back- 
ground process with EPA to reproduce the W(= 
M^) and c.m.-scattering-angle distributions. We 
assume that the angular distribution is the same 
as that in the no-tag process |23| . inspired by the 
observed distribution. The W distribution is tuned 



to the experimental data. Figure [16] shows the ex- 
perimental data and the background MC in a two- 
dimensional plot of the 7r°7r° invariant mass and co- 
sine of the scattering angle (6*) of the 7r° in the 
7r°7r° rest frame with respect to the 77* collision 
axis. The Q 2 dependence of the cross section is tuned 
to reproduce the data before determining the back- 
ground contribution for each Q 2 bin. We normal- 
ize the background size to the observed background 
yield. We find that events from the signal process, 
(e)e7r° with misidentified photons, produce a back- 
ground in the sample of (e)e7r°7r° at extremely for- 
ward or backward scattering angles near an invari- 
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FIG. 13: Results of the fit to the experimental M 77 distributions for representative Q 2 bins. The caption in each 
subfigure shows the type of tag (p- or e-tag), the required trigger, whether or not the Bhabha mask was applied (upper 
line), whether or not the mass and resolution parameters (m and a) are fixed in the fit and the Q 2 range (lower line). 



ant mass of 1 GeV/c 2 . We subtract this component 
in the evaluation of the background. In total, 1305 
events are found for the process e + e~ — > (e)e7r°7r° 
with | cos 0*| < 0.9. 

Using e + e~ — > (e)e7r°7r° MC, the background con- 
tribution from this process is estimated to be about 
2% of the signal yield, with a negligible Q 2 depen- 
dence. This is accounted for as a background factor 



in the calculation of the cross section for the signal 
process. 



4- Background from e + e — > (e)ep°/uj, p°/u) — ► 7r°7 

A similar background from e + e~ — > (e)ep /u>, 
P°/uj —> 7r°7 is also studied. This process can 
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FIG. 14: (a) The -E ra tio distribution for signal candidates for right-sign (RS, dashed line) and wrong-sign (WS, solid 
line) events, (b) The difference between RS events and WS events, (c) The -Bratio distribution for the signal MC 
sample, (d) The M 77 distribution of wrong-sign events with the best fit superimposed (solid line) where the 7r° peak 
and the background contribution (dashed line) are taken into account in the fit. 



contribute to background through a virtual pseudo- 
Compton process, where the final-state photon is re- 
placed by a neutral vector meson. We have searched 
for events that could be described by this process 
by requiring an extra detected photon with a loose 
energy constraint. We find about 180 events with 
the 7r°7 invariant mass peaking at the nominal ui 
mass. Figure [17] shows the experimental data in a 
two-dimensional plot of the 7r°7 invariant mass ver- 
sus the cosine of the scattering angle (8*) of the 7r° 
in the 7r°7 rest frame with respect to the 7r°7 mo- 
mentum direction. Events with an u> dominate, al- 
though it is difficult to separate them from a pos- 
sible p° contribution. In Fig. 17 the horizontal 



band seen at | cos 9* \ > 0.9 is from the signal pro- 
cess (e + e _ — > (e)e7r°) with a misidentified photon. 
This contribution is removed when we estimate the 
p° /u> — > 7r°7 background. 

We simulate this background in MC by replacing 
the e7 system in Rabhat by an ew system with the 
same invariant mass and let them scatter at the same 
angle in the (thus, ew) rest frame. We normal- 
ize the MC sample to the observed number of back- 
ground events and estimate the background to the 
signal process. The result is Q 2 dependent because 
the background process is not a two-photon process; 
it varies from 0.8% for the low Q 2 bin to 3% for 



the highest Q 2 bin. This is accounted for as a back- 
ground factor in the derivation of the cross section 
for the signal process. 



The background yield can be estimated by 
a GVDM (generalized vector-meson dominance 
model) factor in which a photon is replaced by a neu- 
tral vector meson, with branching fractions (BF's) 
for p°/uj — > 7r°7 decays. The GVDM factors for the 
p° and w are assumed to be 1/300 and 1/3000, re- 
spectively [24]. Multiplying them by the BF's to 
77T°, we obtain a factor of 1/32000 for the probabil- 
ity of the background process compared to the VC 
process. With the estimated efficiency of 1.5 x 10~ 4 
from the background MC, the expected background 
yield from the GVDM is 120 events. This estimate 
is within a factor of two of the observed background. 



We conclude that the backgrounds discussed in 
subsections B.l and B.2 are negligible. The total 
contribution from the backgrounds described in B.3 
and B.4 is shown in column 6 of Table |T] 
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FIG. 15: The _E 77 A8 distributions for the (a,d) experimental, (b,e) signal MC, and (c,f) radiative Bhabha background 
MC samples, for two different Q 2 regions indicated in each subfigure. The dashed line shows the selection boundary. 
Arrows show the accepted regions. The background MC sample corresponds to an integrated luminosity 140 times 
smaller than the experimental data. 
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events (a) for the experimental data and (b) for the 
(e)e7T vT° MC samples. 



V. MEASUREMENT OF DIFFERENTIAL 
CROSS SECTION AND TRANSITION FORM 
FACTOR 

In this section, we first describe the method of 
Q 2 unfolding and determination of the representa- 
tive Q 2 value for each Q 2 bin. We then outline how 
the cross section is measured. This is followed by 
the evaluation of the systematic uncertainties and 
some cross-checks. The transition form factor for ir° 
is then determined. Finally, the observed TFF is 
compared with previous experiments and empirical 
parameterizations . 



A. Procedure for Q unfolding 



The 7T° yield is measured in bins of reconstructed 
momentum transfer Q 2 CC . We approximate the true 
value of Q 2 with the corrected value Q 2 OI . We use 
signal MC events to calculate the folding matrix 
for a total of 15 bins in Q 2 , of 1 GeV 2 width for 
5- 12GeV 2 , 2 GeV 2 width for 12-20 GeV 2 , 5 GeV 2 
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FIG. 17: Scatter plot for the selected (e)e7r°7 events: 
| cos 0* | vs. invariant mass of the 7r°7 (see the text for 
the definition of the 9*). 



width for 20 - 30 GeV 2 , and 10 GeV 2 width for 
30 — 50 GeV 2 . This is done separately for the p- 
and e-tags. In our algorithm, we include bins that 
are not used to calculate the TFF. The generated 
events are weighted to have a realistic Q 2 depen- 
dence da/dQ 2 oc Q~ 7 , which is consistent with the 
data. 

The number of events in the same bin before and 
after the unfolding is expressed by 



= > A 



H,rec — / ^ - 
3 



(17) 



where n i lcc is the number of events in the i-th bin 
before the unfolding, which is based on the measured 



Q , and 



''j, cor 



is the number of events in the j-th bin 



after the unfolding, which is based on the true Q . 
By definition, JV Aij is normalized to unity. A^ is 
close to a unit matrix, and the majority of the off- 
diagonal elements are close to zero except for the 
components adjacent to the diagonal. We calculate 
the inverse matrix, which is also close to diagonal. 
The edge bins in the inverse matrix are irregular be- 
cause no experimental measurements are available. 
To avoid bias, we do not derive rij jCOT near the edge 
bins, and use the unfolding method described above 
to obtain results only for the range 7 — 30 GeV 2 . 

For the other bins near the edge, we correct the 
measurement by the factor fa — JV ^ji I Tlj ^ij > 
where Nij is the unnormalized original matrix ob- 
tained from the <5 _7 -dependent MC simulation; this 
method yields correct results when MC events are 
generated with a realistic Q 2 dependence. Figure 18 
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FIG. 18: The number of events for the signal candidates 
in Q 2 bins before and after the unfolding. Note that 
the yield is not normalized to the bin width, and the 
background yield is not subtracted. The open squares 
(dots with error bars) are the yields before (after) the 
Q 2 unfolding. The vertical error bars are statistical and 
are only shown for the plots after the unfolding. The 
data points at Q 2 < 6 GeV 2 for the e-tag are the results 
from the "HiE+50*CsiBB" sample. 



shows the number of events in each bin before and af- 
ter the unfolding. The effect of the unfolding is small, 
but the statistical uncertainty increases slightly. The 
resulted numbers after the unfolding are tabulated as 
V cor in Table m 



B. Representative Q 2 values for each bin 

We define the representative Q 2 for each bin with 
a finite bin width, Q 2 , using the formula 



— (Q 2 ) 
dQ 2 W ' 



1 



AQ 2 



da 
dQ 2 



{Q 2 )dQ\ 



(18) 



where AQ 2 is the bin width. We assume an approx- 
imate dependence of da/dQ 2 oc Q~ 7 for the calcu- 
lation. The representative values are shown in the 
leftmost column of Table HI 
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C. Extraction of differential cross section 

The number of signal events is obtained by the fit 
procedure shown in Sec. IV A in each Q 2 bin, sepa- 
rately for the e-tag and p-tag samples. To take into 
account event migration among bins due to finite 
resolution, the number of events observed in each 
Q 2 bin is corrected using the unfolding method de- 
scribed in Sec. IV Al The number of events is then 
converted to the cross section for e + e~ — > (e)eir° 
using the following formula: 



<$T100 

<D 

CD 



da 
dQ 2 



N(l - b) 



r 1 BfCdt(l + S)AQ 2 



(19) 



where N is the number of signal events in the bin, b is 
the fraction of the estimated ir° peaking background 
in the yield, and rj is the efficiency corrected for 
the beam energy difference between the T(45) and 
T(55 f ), which is the mean of the efficiencies weighted 
by the ratio of the integrated luminosities and the 
2A(Q 2 ) parameter for different beam energies. Us- 
ing this procedure, the differential cross section is ob- 
tained at the e + e~ cm. energy of the T(AS). Here, 
B is the branching fraction for 7r° — > 77, J Cdt is the 
total integrated luminosity, and 1 + 6 is the radiative 
correction for r™ ax — 0.25 used to determine the 
cross section for the tree-level diagrams. The value 
of S is —1% or +2% without or with a correction 
for hadronic vacuum polarization |25j . respectively; 
we adopt the latter value. We include its Q 2 depen- 
dence, which is around 1% in the kinematical region 
of the present measurement, into the systematics. 
The quantity AQ 2 is the bin width in Q 2 . 

We extract the cross section of e + e~ — > (e)e7r° de- 
fined in the range Q\ < 1 GeV 2 , where Q 2 is the 
Q 2 of the photon with smaller virtuality. Our mea- 
surement is performed in a small Q 2 region, mainly 
less than 0.01 GeV 2 , and then the upper bound is 
extrapolated to 1.0 GeV 2 . The requirement for the 
incident photon with the smaller virtuality (which is 
the Q 2 requirement in this paper) depends on the 
experiment and hence the definition of the cross sec- 
tion is somewhat dependent on the experiment. The 
result of this procedure is shown in Fig. [19] The dif- 
ferential cross sections are shown separately for the 
p- and e-tags in Table |TT] 
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FIG. 19: Results for the e+e" -> (e)e-iT differen- 
tial cross section. "Combined" corresponds to twice the 
weighted average of the p-tag and e-tag results. Error 
bars show statistical uncertainties only. The horizontal 
positions in the p-tag and e-tag plots are shifted slightly 
relative to each other to make it easier to see the sizes of 
the statistical uncertainties. 



Fig. [20] We combine the p- and e-tag results for 
the same Q 2 bin, if available, by taking the weighted 
average of the results with a weight that is the in- 
verse square of the statistical error. We then double 
the average. Where only the result from the e-tag 
is available (for the Q 2 range of 4 — 6 GeV 2 ), we 
simply double the value. The differential cross sec- 
tions combined for the p- and e-tags are shown in 
Table HTT1 

The obtained differential cross sections are com- 
pared with those from the previous experiments [TJ U] 
in Fig. [21] Although the definitions of the differen- 
tial cross sections differ slightly between the experi- 
ments, the difference should be compensated in the 
algorithms used to extract the transition form factor. 
Later in this section, we compare in more detail the 
7r°-TFF results obtained by different experiments. 



D. Determination of the transition form factor 

The e + e~ — > (e)e7r° cross section is converted to 
the transition form factor squared by dividing it by 
2A(Q 2 ). Finally, we obtain 



We find that the differential cross sections are con- 
sistent between the p- and e-tag measurements. The 
difference between the two measurements in the same 
Q 2 region is smaller than 1.8tr, where a is the com- 
bined statistical error, for each Q 2 bin as shown in 



Q 2 \F(Q 2 )\ = Q 



2 da/dQ 2 (Q 2 



2A(Q 2 ) ' 



(20) 



where Q 2 replaces the symbol Q 2 . The results for 
this product are shown in Fig. [22] and Table |III| 
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TABLE II: e + e" -> (e)e7r° differential cross section and systematic uncertainties in bins of Q 2 summarize d for the 
p- and e-tags separately. The "Q 2 -indep. syserr" and "Q 2 -dep. eff. syserr" uncertainties are described in Sec 
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Q 2 -dep. eff. 


Fit 


Peaking-bkg 


Total 


Tag (GeV 2 ) 


syserr (%) 


syserr(%) 


syserr (%) 


syserr (%) 


syserr(%) 


6.47 


15.85 ± 1.63 


6 


4 


4 


2 


8 


7.47 


9.47 ±0.76 


6 


2 


4 


2 


8 


8.48 


7.09 ±0.47 


6 


2 


4 


2 


8 


9.48 


4.13 ± 0.29 


6 


2 


4 


2 


8 


10.48 


2.81 ± 0.23 


6 


2 


4 


2 


8 


11.48 


2.31 ± 0.23 


6 


3 


4 


2 


8 


p 12.94 


1.37 ±0.12 


6 


10 


4 


2 


12 


14.95 


0.857 ±0.129 


6 


12 


4 


2 


14 


16.96 


0.539 ±0.127 


6 


10 


4 


2 


12 


18.96 


0.397 ±0.093 


6 


10 


4 


2 


12 


22.29 


0.240 ±0.046 


6 


5 


10 


2 


13 


27.33 


0.172 ±0.042 


6 


5 


10 


2 


13 


34.46 


0.019 ±0.016 


6 


5 


10 


4 


13 


4.46 


37.5 ± 11.2 


6 


2 


7 


2 


10 


5.47 


21.6 ±4.8 


6 


1 


7 


2 


9 


6.47 


15.03 ±2.26 


6 


8 


7 


2 


12 


7.47 


6.45 ± 1.63 


6 


5 


7 


2 


11 


8.48 


6.21 ± 1.02 


6 


5 


7 


2 


11 


9.48 


5.53 ±0.76 


6 


5 


7 


2 


11 


10.48 


2.97 ± 0.48 


6 


5 


7 


2 


11 


e 11.48 


1.77 ±0.49 


6 


5 


7 


2 


11 


12.94 


1.05 ±0.29 


6 


8 


7 


2 


12 


14.95 


0.886 ±0.213 


6 


10 


7 


2 


14 


16.96 


0.607 ±0.183 


6 


10 


7 


2 


14 


18.96 


0.495 ±0.157 


6 


10 


7 


2 


14 


22.29 


0.173 ±0.061 


6 


10 


10 


2 


15 


27.33 


0.088 ±0.051 


6 


10 


10 


2 


15 


34.46 


0.064 ±0.024 


6 


10 


10 


4 


16 



E. Systematic uncertainties and some checks 

The estimates of systematic uncertainties for the 
Q 2 differential cross sections are described in this 
subsection. They can be divided into Q 2 -dependent 
and -independent contributions. 

1. Estimates of systematic uncertainties 

The Q 2 -independent contributions in the system- 
atic uncertainty include tracking uncertainty (1%), 
electron-ID (1%), 77-pair reconstruction efficiency 
(3%), kincmatical selection (2%), geometrical selec- 
tion including the trigger (2%), beam-background ef- 



fects that affect the event reconstruction (2%), and 
integrated luminosity (1.4%). These are correlated 
Q 2 -independent uncertainties, resulting in a total 
of 5% in the quadratic sum for the systematic un- 
certainty arising from the efficiency. Note that we 
do not use a conventional ir° reconstruction from 
a mass-constrained fit. as we fit experimental 77 
mass spectra. The peaking-background contamina- 
tion (l%-4%), fit error in the 7T°-yield extraction 
(5%-10%), and the trigger efficiency (2%-12%) have 
significant Q 2 dependence. We estimate the fit un- 
certainty by taking the difference between the central 
values and a simpler fit with a single Gaussian and a 
linear background applied for the narrower M 77 fit 
range, 0.08 — 0.20 GeV/c 2 . We evaluate the system- 
atic uncertainty due to this M 77 fitting by averaging 
the deviations in the fit results in each of two Q 2 
regions (below and above Q 2 = 20 GeV 2 ) for each 
of the p- and e-tags to compensate for the statistical 
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FIG. 20: Differences between the e + e~ — > (e)e7v° differ- 
ential cross section measured by the p-tag (squares) and 
e-tag (circles) from their weighted mean (equal to a half 
of the combined cross section). The vertical scale is nor- 
malized by the error of the weighted mean (which is the 
same as half of the error for the combined cross section) . 
The horizontal positions of the p-tag and e-tag plots are 
displaced intentionally so that the error bars are visible. 



> 
CD 

g 100 

CM 

a 
~& 

73 10 



0.01 



% 



X 



• Belle 

♦ CLEO 
BaBar 



10 



20 



30 

Q 2 (GeV 2 ) 



40 



FIG. 21: Comparison of the results for the e + e~ — > 
(e)en° differential cross sections among different experi- 
ments. The vertical error bars are statistical only. The 
algorithms used to measure the cross section differ be- 
tween the experiments. 



fluctuations, as listed in column 6 of Table |TT] 

To estimate the systematic uncertainty in the trig- 
ger efficiency, we study the efficiency using a radia- 
tive Bhabha process; we find agreement between the 
data and the Rabhat MC at the 5 — 10% level, as 

In addition, 



shown in Fig. 30'c,d) (Appendix B3). 
the comparison of the absolute yields of the radia- 
tive Bhabha events shows consistency between the 



TABLE III: e+e" -> (e)e7T° differential cross section 
combined for the p- and e-tags with systematic uncer- 
tainties (e sy s) and the transition form factor Q 2 \F(Q 2 )\ . 
The first and second uncertainties for Q 2 \F(Q 2 )\ are sta- 
tistical and systematic, respectively. 
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FIG. 22: Results for the tv° transition form factor from 
the present measurement multiplied by Q 2 . The error 
bars are statistical only. The dashed line shows the 
asymptotic limit from pQCD (~ 0.185 GeV). 



data and MC also within 10% level, which is com- 
parable to the estimated size of systematic uncer- 
tainty without including the uncertainty from the 
trigger efficiency of ~ 8% (see Fig. [29|c,d) and the 
text in Appendix B 3 1 . Therefore, we conclude that 



the systematic uncertainty in the trigger efficiency 
does not exceed the deviations seen in studies us- 
ing radiative Bhabha events: 12% for the p-tag and 
10% for the e-tag, where the values are taken from 
the overall tendency of the deviations in the Q 2 de- 
pendence shown in Figs. |29|c,d) and(307c,d). How- 
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ever, we do expect the systematic uncertainty to be 
smaller for many individual Q 2 bins, e.g., in bins 
where the effect of the Bhabha-veto threshold is rel- 
atively small and the trigger efficiency is not very 
sensitive to it. To obtain the Q 2 -dependent system- 
atic error, we have varied the Bhabha-veto trigger 
threshold by ±0.2 GeV in the trigger simulation ap- 
plied to the signal MC events and taken the change 
of the efficiency for each bin if it is smaller than the 
maximum values mentioned above. In practice, we 
have found that the variation has a large Q 2 depen- 
dence, ranging from 2% to 25%, in each of the p- 
and e-tag samples. We replace the uncertainty deter- 
mined by the variation of the trigger threshold in the 
signal MC by 12% (10%) for the p-tag (e-tag) case if 
the former is greater than 12% (10%) in an individ- 
ual Q 2 bin. We replace the systematic uncertainty 
in the bin Q 2 = 14 — 16 GeV 2 for the p-tag sample 
and in the six bins in the Q 2 range 14 — 40 GeV 2 for 
the e-tag sample. 

We assign a ±(1% — 3%) systematic uncertainty 
from a possible deviation of the Q 2 dependence in the 
signal MC from the data distribution for the lowest 
two bin results in each of the p- and e-tags. In other 
Q 2 regions, this uncertainty is negligibly small. We 
take the difference between the calculated efficiencies 
with two assumed Q 2 dependences: |i^(Q 2 )| 2 ~ 1 
and ^(Q 2 )! 2 ~ 1/Q 2 - The combined uncertainties 
from this source and from the trigger emciency are 
shown in Table |ll| as a Q 2 -dependent systematic un- 
certainty in the emciency. The evaluation of uncer- 
tainties is performed separately for the p- and e-tags. 

In addition, the uncertainties from the calculation 
of the radiative corrections (3%) and assumed form 
factor effects for the photon with smaller virtuality 
(1.0%) enter when determining the differential cross 
section. The correlated Q 2 -independent uncertainty 
is in total 6%, of which 5% is from the efficiency de- 
termination as described above, in a quadrature sum 
for the uncertainty in the differential cross section. 

The total systematic uncertainties for the com- 
bined differential cross section and for the transi- 
tion form factor are shown in Tables |TT] and IIIII We 
assume that the systematic uncertainties in the p- 
and e-tags correlate maximally. Because the tran- 
sition form factor is derived from the square root 
of the cross section, the relative systematic uncer- 
tainty is half of the above. The calculation uncer- 
tainty for the conversion factor 2A(Q 2 ), 2%, is small 
compared to other sources of uncertainties. The un- 
certainty in the value of Q 2 in each bin is also neg- 
ligible (AQ 2 /Q 2 ~ 0.3%) because Q 2 \F{Q 2 )\ varies 
slowly with Q 2 . 



2. Check for ISR in the signal process 

We have included the effect of ISR in the signal 
MC simulation but only for the tagged electron side 
with exponentiation. This is a reasonable approxi- 
mation because only the tag side shifts Q 2 and kine- 
matics significantly, while low-Q 2 radiation from the 
untagged side does not change any observable. The 
corrections in S are small for both cases (at the 2% 
level) . We compare the r k distributions in data after 
background subtraction using the 7r°-mass sideband 
events and the signal MC samples after the selection 
as shown in Fig. [23j The fraction of energy from 
radiation relative to the beam energy, (denned in 
Eq.(14l), is calculated by an approximate kinemati- 
cal relation from the observed variables, 

r k = ^-(K + E;)-\K+Pt\. (21) 

In the BaBar analysis, this parameter is directly used 
as a selection pQ . In our analysis, instead of ru we use 
the variable i?ratio- These two variables are strongly 
correlated, as, kinematically, -E ra tio = 1 when = 0. 
The peak and the core part of the r k distribution 
(—0.02 < rfe < 0.06) are safely contained by our se- 
lection criteria for -E ra ti OI so the consistency between 
the data and MC after this selection permits a re- 
liable determination of the signal efficiency. As the 
signal MC events are generated with r™ ax = 0.25, no 
bias is introduced in the signal region by our selec- 
tion criteria. 



F. Comparison with previous measurements 
and empirical fits 

The measured product Q 2 \F(Q 2 )\ is compared 
with the previous measurements [U |3j 0] in Fig. |24| 
Our results agree with previous measurements within 
uncertainties for Q 2 below 9 GeV 2 ; in the higher Q 2 
region, they do not show the rapid growth observed 
by BaBar pQ. 

BaBar has suggested a parameterization fit with 
the functional form 



Q Z \F(Q')\=A[ y „ 
1 10 GeV 



(22) 



where A and (3 are fit parameters. BaBar reported 
A = 0.182 ± 0.002 GeV and /3 = 0.25 ± 0.02; how- 
ever, the uncertainty for A does not include the Q 2 - 
independent systematic uncertainty that amounts to 
2.3% (namely ±0.0041 GeV as this error component 
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FIG. 23: Distribution of the ISR energy fraction (r^) 
calculated from observables in data (solid dots) and in 
the signal MC samples (histogram). The data distribu- 
tion is shown after background subtraction. The MC 
events are normalized to the experimental yield after 
background subtraction. Note that there are uncertain- 
ties in the value of rk originating from the measurement 
and kinematical approximation. 
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FIG. 24: Comparison of the results for the product 
Q 2 \F(Q 2 )\ for the n° from different experiments. The 
error bars are a quadratic sum of statistical and system- 
atic uncertainties. For the Belle and BaBar results, only 
a Q 2 -dependent systematic-error component is included. 
The two curves denoted fit(A) use the BaBar parameter- 
ization while the curve denoted fit(B) uses Eq.( |23[ ) (see 
the text). The dashed line shows the asymptotic predic- 
tion from pQCD (~ 0.185 GeV). 



for Q 2 \F(Q 2 )\). Taking this into account, the uncer- 
tainty on A reported by BaBar could be replaced by 
A = 0.182 ± 0.005 GeV. To compare our results with 
BaBar's, we use the same parameterization in our fit 
procedure and assume a Q 2 -independent systematic 
uncertainty (thus, the total normalization error) of 
3.2% for Q 2 \F(Q 2 )\; we remove this component from 
the combined statistical and systematic errors and 
instead add it in quadrature to the uncertainty in A. 
The fit results from Belle are A = 0.169 ±0.006 GeV 
and /3 = 0.18 ± 0.05. The goodness of the fit is 
X 2 / n df = 6.90/13, where ndf is the number of de- 
grees of freedom. The fit results are also shown in 
Fig. [24] The fit of the Belle data to the function is 
good, and we find a difference of ~ 1.5cr between the 
Belle and BaBar results in both A and /?. 

We then try another parameterization in which 
Q 2 \F(Q 2 )\ approaches an asymptotic value, namely 

q2|f(q2)I = o?tV (23) 

The fit gives B = 0.209 ± 0.016 GeV and C = 2.2 ± 
0.8 GeV 2 with x 2 /ndf = 7.07/13, and is also shown 
in Fig. [24] The fitted asymptotic value, B, is slightly 
larger than the pQCD value of ~ 0.185 GeV but is 
consistent. 

For a simple estimate of the consistency between 
the Belle and BaBar results, we compare the data 
from individual experiments with a reference curve 



obtained by fitting the data from both experiments 
together using parameterization of Eq. (23). The 
seven data points from Belle (from BaBar) for the 
range 9 GeV 2 < Q 2 < 20 GeV 2 , where the two mea- 
surements seem to be systematically shifted, devi- 
ate from the reference curve by (—6.1 ± 3.8)% (by 
(+4.8 ± 3.0)%) in average relative to it. We incor- 
porate the Q 2 -independent uncertainty in each mea- 
surement in the above error. The difference between 
the Belle and BaBar deviations, (10.9 ± 4.8)%, cor- 
responds to a 2.3cr significance. This result does not 
depend on choice of the reference curve. 



VI. SUMMARY AND CONCLUSION 

We have presented a measurement of the neutral 
pion transition form factor for the process 77* — > 7r° 
in the region 4 GeV 2 ~ Q 2 ~ 40 GeV 2 with a 
759 fb _1 data sample collected with the Belle detec- 
tor at the KEKB asymmetric-energy e + e~ collider. 
The measured values of Q 2 \F(Q 2 )\ agree with the 

previous measurements [TJ 02 0] for Q 2 ^9 GeV 2 . 
In the higher Q 2 region, in contrast to BaBar, our 
results do not show a rapid growth with Q 2 and are 
closer to theoretical expectations [5]. 
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Appendix A: Bhabha-veto and Bhabha-mask 
selection 

The Bhabha veto in the HiE trigger of the Belle 
trigger system rejects a significant portion of the sig- 
nal |26| . An accurate determination of the efficiency 



from this veto is essential for a reliable measurement 
of the signal cross section. In this section we provide 
a detailed description of the Bhabha-mask rejection 
that selection criterion (9) is based on. 



1. Bhabha-veto trigger 

A detailed description of the Bhabha trigger 
(CsiBB) at Belle is given elsewhere [TS]. It fires 
when the energy sum over one- or two-"#-ring" 
groups among eleven possible combinations (referred 
to hereafter as "Logic") is higher than a given thresh- 
old. A 0-ring is an array of CsI(Tl) crystals in a 
certain polar angular range covering 2n in azimuth. 
A #-ring combination corresponds to back-to-back 
Bhabha electrons in that angular range. 

The sum of energies over two 0-rings is required 
but the energy can concentrate in either ring. Thus, 
an event with a single cluster can satisfy the trigger 
logic even if it has a non-back-to-back topology in 
the e + e~ cm. frame. This vetoes some fraction of 
the signal in the two-dimensional angular acceptance 
plane for (cos# e , cos# 77 ), especially in the forward 
region where the energies of the electron and photons 
are relatively high. 



2. Bhabha-mask rejection 

As mentioned in selection criterion (9), we reject 
events for which the angular combination (cos(9 e , 
cos# 77 ) falls in the predetermined two-dimensional 
angular area. We divide the polar angle region 
for the selection of an electron and photon system, 
-0.6235 < cos6» < 0.9481, that is, 18.53° < 9 < 
128.57°, into 14 subregions with 13 boundaries [57] 
according to the boundaries of the Csl crystals that 
belong to different trigger segments. The acceptance 
plane (cos6* e , cos(9 77 ) is then subdivided into 14 x 14 
rectangular cells (see Fig. [2j. 

We estimate the efficiency of the Bhabha veto in 
each cell using the trigger simulator, TSIM. The cells 
that are sensitive to the veto are removed from the 
selection. The accepted cells are shown in Fig. [2] 
in yellow. In this figure, we overlap the distribu- 
tion from the signal MC simulation (shown with 
the scattered dots), where the Bhabha veto is not 
applied, but the high-Q 2 contribution is enhanced 
as compared to the actual distribution. The two- 
dimensional angular regions for the p-tag and e-tag 
are separated into the left-lower region and the right- 
upper region, respectively. 
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3. Unbiased sample 

The HiE trigger is vetoed by the CsiBB trigger. 
Hence, the "HiE OR CsiBB" data sample, which is 
triggered by either of them, is completely free from 
the Bhabha-veto logic. The events triggered by the 
Bhabha (CsiBB) trigger are recorded after prescaling 
by a factor of 1/50. These prescaled events are used 
in the analysis of the ir° TFF in the region of Q 2 
between 4 and 6 GeV 2 . 

In addition, this sample is very useful to calibrate 
the energy threshold of the CsiBB trigger as well 
as to test the trigger simulation of the HiE trigger 
in which the Bhabha veto is included. These cali- 
brations and the test of the trigger simulation have 
been carried out by using radiative Bhabha events 
(e + e~ — > (e)e"/) in the special configuration where 
one of the final state electrons is scattered into the 
beam pipe. This configuration is usually referred 
to as Virtual Compton scattering (VC). The kine- 
matics of VC are very similar to the signal process 
(e + e~ — > (e)e7r°) so we can use this sample as a cali- 
bration signal to validate the trigger simulation using 
data, which is one of the most important aspects of 
this analysis. 

To compensate for prescaling, we multiply by 50 
the number of events recorded by the CsiBB trigger. 
Thus, to emphasize that the effect of the Bhabha 
veto is statistically compensated in the complete 
dataset, it is referred to as "HiE+50*CsiBB". 



Appendix B: Calibration of the trigger threshold 

It is necessary to calibrate the trigger thresholds 
for the Bhabha veto because the inefficiency due to 
the trigger for the signal is significant. In contrast, 
the threshold for the HiE trigger is efficient because 
it is low (around 1.15 GeV) compared to our analysis 
selection requirement . 

There are eleven trigger thresholds for the Bhabha 
veto (the same as the CsiBB trigger) and their nom- 
inal value is 5.0 GeV, except for Logic#2 (5.5 GeV) 
and Logic#ll (3.0 GeV) [TS]. However, the thresh- 
olds determined from Bhabha events in data indicate 
that the actual values of the thresholds are larger, 
typically by 0.2 GeV. In addition, the threshold im- 
plemented in the TSIM trigger simulator is a step 
function while the data exhibit the turn-on curves 
typical for such triggers. We correct for these differ- 
ences between data and MC. 

The trigger efficiency is studied and tuned using 
radiative Bhabha (e + e~ — > (e)e^f) events with the 
VC configuration as described in the next subsec- 



tion. The radiative Bhabha cross section in the VC 
configuration is very large; it is O(10 4 — 10 5 ) times 
larger than the signal process. Significant and suffi- 
cient statistics for radiative Bhabha events are avail- 
able even in the prescaled CsiBB data samples. 

In Sec. |B 1[ we describe tuning of the Bhabha- 
veto threshold using an unbiased sample of radiative 
Bhabha events, which are compared to events gener- 
ated by Rabhat MC together with TSIM and GSIM 
(detector simulator). The Rabhat MC is described 
in detail in Sec. |B 2| Finally, in Sec. |B 3| various 
distributions including Q 2 in data and MC are com- 
pared to check the validity of the tuned MC, first for 
unbiased data and then for HiE triggered data. 



Tuning of CsiBB thresholds with radiative 
Bhabha events 



We compare the radiative Bhabha events in data 
with the MC events from Rabhat, collected with se- 
lection criteria similar to those for the signal pro- 
cess. To select the radiative Bhabha events, we relax 
the selection requirements by allowing events with 
a single photon. The energy requirement for both 
the electron and photon is slightly more restrictive, 
greater than 1.5 GeV, in order to collect (e)e7 sam- 
ples with negligible background contamination. The 
conditions for the angles and the other selection re- 
quirements are the same as those in the analysis 
of the signal process. The three-body kinematical 
selection is applied assuming a single real photon, 
0.85 < £ratio < 1-1. We define the p- and e-tags and 
Q 2 in the same manner as in the measurement of the 
single-tag two-photon process from the detected elec- 
tron for each event, although the physics mechanism 
is different . 

We measure the threshold curve for each CsiBB 
Logic by counting the yields from the two triggers 
and by forming the ratio 



2/ratio 



V 



50 * jVpsiBB 

50 * N, 



HiE 



CsiBB 



(Bl) 



as a function of the energy sum (-E SU m) over the elec- 
tron and photon deposited in a pattern of the Logic 
for both experimental data and MC, where ./VcsiBB 
(-VHiE+50*iVc s iBB) is the number of events triggered 
by the CsiBB (HiE + CsiBB multiplied by 50) trig- 
ger^) and the factor of 50 corrects for the prescaling 
in the CsiBB trigger. Note that, in matching ?/ ra tio 
between data and MC, tuning of the TSIM threshold 
is essential while the details of the physics simulated 
by Rabhat are not. We obtain the energy depen- 
dence of y r atio for each Logic. We use the energy 
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FIG. 25: The threshold behavior for the four selected Logics of the CsiBB trigger (identified in each subfigure) obtained 
from the experimental samples (left column) of radiative Bhabha events and from the MC events, where the trigger 
simulator is tuned (right column). The dashed lines show the position of the nominal trigger thresholds. 



sum (-Esum) over the electron or the photon that en- 
ters the relevant ECL #-ring(s) to the Logic. In each 
Logic, only the events in which the sum of the de- 
posited energy to the corresponding Logic exceeds 
•E-max — 1 GeV are used, where E max is the greatest 
energy deposit in a given Logic (among the eleven) 
in each event. 

We obtain the threshold curves for each Logic for 
data and MC simulation and tune the latter to match 
the former. We tune the threshold behavior of the 
MC so that it agrees well with that in data as shown 
in Fig. [25] The thresholds are thus adjusted with 
an accuracy better than ±0.2 GeV. The non-zero ef- 
ficiency below the nominal threshold is due to an 
artifact from events that are triggered by a differ- 
ent Logic. The real thresholds are near 5-6 GeV for 
Logics#l — 10 and near 3 GeV for Logic#ll. The 
gaps near 7.5 GeV, when present, correspond to a 
region almost prohibited by the kinematics of (e)e7- 

We confirm the result of the tuning and its validity 
by comparing distributions of the radiative Bhabha 
events between the experimental data and MC in 
Sec. |B 3| We use the event yields for this purpose 



because the radiative Bhabha cross section is well 
known. 



2. Rabhat, MC code for radiative Bhabha 
events 

We use the MC simulation program Rabhat [T6"| 
to generate events with a virtual Compton configu- 
ration. Rabhat simulates the radiative Bhabha pro- 
cess with a i-channel mass singularity. This singu- 
larity occurs in a topology in which the virtual pho- 
ton that is emitted at an extremely forward angle 
from one of the incident electrons is hard scattered 
off the other electron. Although the cross section 
is finite and is not very large, our requirement of 
an energetic final-state photon at a finite angle re- 
sults in a relatively sharp mass singularity for an 
acoplanarity angle of order E 1 /m e . Rabhat adopts 
a careful choice of integration variables and a proper 
technique for numerical integration. We used Rab- 
hat to simulate approximately 2 x 10 s (3.6 x 10 7 ) MC 
events for yf~s = 10.58 GeV (10.88 GeV). A limita- 
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FIG. 26: Comparison of the experimental radiative- 
Bhabha samples (dots) and the corresponding MC (Rab- 
hat, histogram) normalized to the integrated luminosity, 
triggered by HiE or CsiBB, where the Bhabha mask is 
not applied. The distributions are (a) polar angle of the 
photon, (b) polar angle of the electron, (c) sum of the 
photon and electron energies and (d) identification num- 
ber of the cell in the (cos# e , cos# 7 ) plane. 



FIG. 27: (a,b) Comparison of the experimental radiative- 
Bhabha samples (dots) and the corresponding MC (Rab- 
hat, histogram) normalized to the integrated luminosity, 
triggered by HiE or CsiBB, where the Bhabha mask is not 
applied, for the Q 2 distributions of the detected electron, 
where the two selection criteria related to the p t balance 
are loosened (see the text). (c,d) The ratio of the exper- 
imental yield to the normalized MC yield. The sizes of 
the vertical error bars correspond to the statistical errors 
in both experimental and MC samples. The horizontal 
lines at 1.0 show the Q 2 region where the measurement 
of the 7T° TFF is performed. 



tion of Rabhat is that it calculates the lowest order 
e+e~ — > e+e~7 process (C(a 3 )). There could be 
an uncertainty from radiative corrections, typically 
of order 5%, which is estimated from the size of the 
radiative correction 1 + S for the Compton process, 
67 — > &y, where 5 ~ —6% in the similar ej cm. 
energy region (around 3-10 GeV) [25] . 



3. Comparison between the data and MC for 
radiative Bhabha events 



In this subsection, we first compare experimental 
data for the HiE+50*CsiBB sample with MC, where 
the CsiBB threshold effect is absent, to show that 
MC can reproduce various distributions fairly well 
(Figs. 



26 and 27 1. Similar plots are presented for 



the HiE triggered events (Figs. K28^ and 29 1 to show 
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FIG. 28: Comparison of the experimental radiative- 
Bhabha samples (dots) and the corresponding MC (Rab- 
hat, histogram) normalized to the integrated luminosity, 
triggered by HiE, where the Bhabha-mask and -veto are 
applied. See the caption of Fig. [26] for more information. 



that the agreement between data and MC is also 
adequate. We use an MC sample where the CsiBB 
threshold is tuned as described in Sec. IB II We then 
show the HiE trigger efficiency determined by data 
using radiative Bhabha events and compare it to the 
results from the trigger simulation (Fig. 30 1. 

We show the comparison in distributions for polar 
angles of the photon and electron, energy-sum and 
cell pattern as correlation of the angles in Fig. [26] 
for the sample of HiE+50*CsiBB. The MC events 
are normalized to the nominal luminosity; radiative 
corrections are included in this normalization. The 
horizontal axis for the cell identification number cor- 
responds to each cell in the (cos6* e , cos(9 7 ) plane 
in which the symmetry in interchanging cos 6 e and 
cos6* 7 is taken into account (see Fig. [2]). 



FIG. 29: Comparison of the experimental radiative- 
Bhabha samples (dots) and the corresponding MC (Rab- 
hat, histogram) normalized to the integrated luminosity, 
triggered by HiE, where the Bhabha-mask and -veto are 
applied. See the caption of Fig. [27] for each figure. The 
horizontal lines at 1.0 in (c,d) show the Q 2 region where 
the measurement of the n° TFF is performed using HiE- 
trigger data. 



The Q 2 distributions for data and MC for the 
unbiased sample (HiE+50*CsiBB) are shown in 
and 



Figs. 27 ^a) and (b) for the p- and e-tag samples, 
respectively. Figures |27{c) and (d) are the ratio be- 
tween data and MC for the p- and e-tag separately. 
In these comparisons, we include the radiative cor- 
rection of S = —6% and loosen the selection criteria 
related to the p t balance conditions for the ej system 
to a(e,7) < 0.3 radians and |Sp^| < 0.5 GeV/c, for 
a more meaningful comparison of the yield, noting 
that the p t -balance is affected by the emission of an 
extra photon; the effect is typically 4%. 

Figures [28] and [29] present similar distributions 
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FIG. 30: (a,b) Comparison between the efficiency of the "Bhabha-mask x Bhabha-veto" for the radiative-Bhabha 
samples (dots, a; rat i (data) described in the text) and the corresponding MC (Rabhat, histogram, i rst io(MC)) for (a) 
p-tag and (b) e-tag events. (c,d) Ratio (sdr) of the efficiencies which are plotted in (a,b) for the radiative-Bhabha 
samples and the corresponding MC for (c) p-tag and (d) e-tag events. The vertical error size is derived from the 
statistical errors in both experimental and MC samples. The horizontal lines at 1.0 show the Q 2 region where the 
measurement of the ty° TFF is performed using the Bhabha mask and Bhabha veto. 



for the samples after applying TSIM for the HiE 
trigger including applications of the Bhabha veto 
and Bhabha mask. The results agree to within the 
5 - 10% level (Fig. |29Tc), (d)). A deficit of events 



is seen in both the p-tag and e-tag modes, and also 
in the HiE+50*CsiBB sample (Fig.^Tjc), (d)). This 
discrepancy is not due to the Bhabha veto, because 
a deficit of similar size is also seen in the regions 
where the veto has little or no inefficiency, e.g., for 
Q 2 < 10 GeV 2 with the p-tag. Moreover, the deficits 
are of similar size in the unbiased sample and in the 
HiE-trigger sample. These deficits could be due to 
the uncertainty of the radiative correction S (4%) , an 
uncertainty in the efficiency due to additional multi- 
photon emission, which is not accounted for in the 
Rabhat MC (4%), as well as the experimental mea- 
surement itself (~ 5%, similar to the Q 2 -independent 
systematic uncertainty in the efficiency for the sig- 
nal process in Sec. 



total. 

The ratio z ra ti , 



V E I , which amounts to 



defined as 



HiE 



50 * N G 



siBB 



7c in 



(B2) 



i.e., the number of events triggered by HiE divided 
by those triggered by "HiE OR CsiBB", provides the 
fraction of the events that survive the Bhabha veto. 
The Bhabha mask is (is not) applied to the HiE 
(HiE+50*CsiBB) sample. Thus, x rat i contains the 
effect of the Bhabha mask as an inefficiency. Note 



that this ratio is not affected by uncertainties in 
normalization such as those coming from the radia- 
tive corrections of the Bhabha generators as well as 
that from the experimental measurements. The ratio 
a^ratio is shown as a function of Q 2 for the data and 
Rabhat MC events tagged by positrons and electrons 
in Fig. 30 ^a) and (b), respectively. The double ratio 
(cdr) between data and MC 



Zratio(data) 
,(MC) ' 



(B3) 



shown in Fig. |30[c) and (d), demonstrates that the 
Bhabha-veto efficiency simulated by the MC agrees 
with data to within ±(5% - 10%). 
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